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Abstract —This paper presents a method developed 
for finding sinusoidal components within a nonlinear 
non-stationary time-series data using Genetic Algo¬ 
rithm (GA) (a global optimization technique). It is 
called Search-Enhanced Instantaneous Frequency De¬ 
tection (SEIFD) algorithm. The GA adaptively define 
the configuration of the components by simulating the 
solution finding process as a series of genetic evolutions. 
The start-time, end-time, frequency, and phase of each 
of these components are identified once convergence in 
the implementation is achieved. 

Index Terms —Adaptive data analysis; Instantaneous 
frequency; Global optimizations; Genetic algorithms. 


consituting sinusoidal components in order to reveal the 
instantaneous frequencies of the time-series data. 

II. The implementation 

The sinusoidal parameters are encoded as a string of 
binary values. The length of the binary string is deter¬ 
mined based on the sum of the bit counts carried by each 
of the parameters computed with their respective lower 
and upper bound, as well as the step size of increment. 
The bit count for each of the parameters is computed as 
below. 


I. Introduction 
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The Search-Enhanced Instantaneous Frequency Detec¬ 
tion (SEIFD) algorithm uses Genetic Algorithm (GA) [3]i 
|5j to identify the characteristics of sinusoidal components 
(see Fig. Q within a nonlinear non-stationary time series 
data. The time-series data is modeled as a superposition of 
these components, each exist within a unique time range. 
It is represented as below. 

N 

f[t) = X Sin ( 2 Kfnt + 2 (f> n t) (1) 
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where, n is the index of each sinusoidal component; a n is 
its amplitude; /„, is its frequency; and cf> n is its phase. 


where, bit(:r) is the bit count of parameter x\ x u b and 
x\b respectively are the upper bound and lower bound of 
parameter x ; and Ax is the step size of each increment 
of parameter x in the solution. The solution string is 
organized as shown in Fig. [2} 
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Fig. 2. A solution string in GA. 

The real value of each of the parameters can be decoded 
from the solution string using the following equation. 



Fig. 1. A sinusoidal component within a time-series data. 
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where, l is the length of bits in the solution string repre¬ 
senting the value of x. 

The solution string is regarded as a chromosome (Fig. 
[d]) in the simulated genetic evolution process. In the initial 
stage, GA generates a group (a population) of random 
chromosomes of varying quality computed with the fol¬ 
lowing fitness function. 


The parameters which define each of the components 
are encoded as a solution string consisting binary values. 
During implementation, the GA attempts to find their 
value by simulating it as a solution improvement process 
modeled as a series of genetic evolutions. 

Because the problem is non-stationary, each of the 
sinusoidal components possesses unique start and end 
time (denoted respectively as t" t and t”). They are as 
well encoded in the solution string. In this way, the 
solver attempts to decompose the time-series data into its 
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where the R{t) is the raw time-series data. 

The fitness function describes the environment (the 
solution space) in which the chromosomes evolve. The 
implementation should achive continuous improvement 
(lower fitness function in every subsequent generation) 
until the fitness score converged at a targeted fitness score. 
The targeted fitness score of F(x) is 0. 
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Fig. 3. The iteration-best value of a n , / n , and 4> n found by the GA. 


Gene The value is called allele 

0 110 1 l^P^O 0010100110110 

Chromosome 

Fig. 4. The structure of chromosome in GA 

Then, GA selects some chromosomes possessing good 
fitness score for creating new offsprings in the subsequent 
repetitions. In the process, the algorithm identify and 
retain good alleles that lead to the yielding of better 
fitness scores. Meanwhile, it modifies the inferior allele 
values using simulated biological crossover and mutation 
processes. The operation is repeated until the termination 
criterion is finally met. 

Compare to Fourier transforms, the method adaptively 
identifies the start and end of such sinusoidal components 
(within the time horizon of the data), as well as their exact 
instantaneous attributes. Among others, its performance 
is not constrained by: i.) the need to ensure the completion 
of at least one oscillation cycle in all of the sinusoidal com¬ 
ponents; and ii.) the existence of discontinuity within the 
data. The result fulfills the following generalized Fourier 
representation jl]: 

X(t) = y^qj(f)exp J coi(t)dt^j (5) 

The identification of the quantity of sinusoidal compo¬ 
nents at every time step, nt can be solved as a Combinato¬ 
rial Optimization (CO) problem (the GA is a well-known 
approach for solving CO problems (i.e., n, H)). To do so, 
the method could be implemented as below: 


1) Establish a list of repeatable combination, C of 
length L, which consist of a list of integer, n t < N*-, 
where L is the length of the data, N* is an integer 
representing the maximum quantity of sinusoidal 
components in the data. 

2) Implement GA to determine the value of components 
in C, Ci. 

3) Determine the t st j and t e l for each of the sinusoidal 
component based on C, and identify their configura¬ 
tions. 

4) Improve C until convergence is achieved. 

III. Test case 

The method is tested on a time series data shown 
in Fig. [5] The algorithm is codified using the Python 
programming language. 



Fig. 5. A sample nonlinear non-stationary time-series data. 

The iteration-best value of a n , f n , and 4> n found by the 
GA is shown in Fig.[3| The changes in the parameter values 
which leads to the convergence of F(x) is clearly seen in 
the figure. The final results obtained in the implementa¬ 
tion is shown in Fig. [6] 
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Fig. 6. The sinusoidal components in the test data. 
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